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Adsorption of molecular oxygen on B-, N-, A1-, Si-, P-, Cr- and Mn-doped graphene is theoretically 
studied using density functional theory in order to clarify if O2 can change the possibility of using 
doped graphene for gas sensors, electronic and spintronic devices. O2 is physisorbed on B-, and N- 
doped graphene with small adsorption energy and long distance from the graphene plane, indicating 
the oxidation will not happen; chemisorption is observed on A1-, Si-, P-, Cr- and Mn-doped graphene. 
The local curvature caused by the large bond length of X-C (X represents the dopants) relative to C- 
C bond plays a very important role in this chemisorption. The chemisorption of O2 induces dramatic 
changes of electronic structures and localized spin polarization of doped graphene, and in particular, 
chemisorption of O2 on Cr-doped graphene is antiferromagnetic. The analysis of electronic density 
of states shows the contribution of the hybridization between O and dopants is mainly from the p 
or d orbit als. Furthermore, spin density shows that the magnetization locates mainly around the 
doped atoms, which may be responsible for the Kondo effect. These special properties supply a good 
choice to control the electronic properties and spin polarization in the field of graphene engineering. 

PACS numbers: 61.48.De, 68.43.-h, 75.20.Hr, 73.22.-f 



I. INTRODUCTION 

Graphene has been becoming a new star in a lot of 
fields after its successful fabrication^^, especially for the 
application of two aspects below because of the many 
reasons for the renewed interest: first, it is very potential 
to apply graphene as gas sensors with high sensitivity 
because the transport properties exhibit large changes 
upon exposure to several gases such as N02^. Graphene 
can be used as an excellent sensor material because of 
its special properties such as two dimensional property 
maximizing the interaction of adsorbates on the layer, 
low Johnson noise^'^^, and few crystal defects^'^'^. This 
tells us the electronic properties of graphene system can 
be sensitive to the adsorption of gases. Second, there 
is an extraordinary interest in the electronic band prop- 
erties and magnetic order in graphene-related materials 
that can be put into use as the next-generation electronic 
devices and recording media, magnetic inks, spin qubits 
and spintronic devices^'^'^, which can use the advantages 
of high mobility of electrons in graphene and long coher- 
ence times in carbon-based materials. A few investiga- 
tions have been performed until now to search for the 
potential application of graphene or graphene nanorib- 
bons as gas sensors^^'^^. However, it has been shown 
that intrinsic graphene can only physisorb most of gas 
molecules^ and has no band gap and spin polariza- 
tion. This, for example, prevents the use of graphene in 
making gas sensors, transistors and spintronic devices. 

In order to improve the sensitivity for gases and the 
electronic structures in graphene and carbon nanotubes 
(CNTs), the method of doping is usually used. B- and 
N-doped graphene can also improve gas sensing of CO, 
NO, NO2 and NHa^^; Al-doped graphene can be sen- 
sitive to most of common gases in air^^; H2O on Ti- 
doped graphene^^, chlorophenols on Si-doped CNTs^^ 



are also reported. Meanwhile, the electronic struc- 
tures and magnetic properties in graphene and graphene 
nanoribbons can be changed by doping different atoms 
such as boron^^, sulfur and phosphorus^^, transition 
metal atoms^^^^, or by defect engineering^^, function- 
alized with different atoms, molecules and clusters^^ and 
so on. Experimentally, several CNTs and graphene- 
based materials with different dopants such as nitrogen^^, 
boron^^, phosphorus^^ were also synthesized. It has been 
shown by these studies that the doped graphene or CNTs 
materials are very potential for gas sensors, electronics 
and spintronics. 

However, in order to use these materials in reality, 
the effect of O2 should be considered, since O2 is one 
of the most important gases taking up more than 20% 
in air. Meanwhile, it has been shown that O2 can hole- 
dope semiconductors^^. Also, gas molecules can change 
the electronic and magnetic properties of graphene and 
graphene nanoribbons^^. Previous calculations^^ and 
experiments^^ have shown that O2 in the triplet state 
is physisorbed on the surface of single-walled carbon 
nanotubes (SWCNTs) and graphene, and the transport 
properties could not be significantly affected. On the 
other hand, the electronic properties of carbon-doped 
boron nitride nanotubes can be changed dramatically by 
the chemisorption of O2 molecule^^. This indicates that 
the dopants can improve the reactivity of materials, and 
the O2 molecules can give rise to very different proper- 
ties. Therefore, how O2 molecule can affect the proper- 
ties of doped graphene is necessary to be understood in 
the fields of both gas sensing and electronics. 

In this paper we study from first principles the adsorp- 
tion of molecular oxygen on B-, N-, A1-, Si-, P-, Cr- and 
Mn-doped graphene. In order to illuminate easily, these 
materials are abbreviated to be EG, NG, AG, SG, PC, 
CG and MG. BG and NG retain a planar form, while 



other atoms protrude out of the graphene layer and in- 
duce a local curvature in graphene. Furthermore, O2 
molecule is physisorbed on BG and NG with relatively 
small adsorption energy and large distance of X-0 (X rep- 
resents the doping atom), comparing with chemisorption 
on the doped graphene with other dopants. In the end, 
the chemisorption induces dramatic change for the elec- 
tronic structures of doped graphene, and injects magnetic 
moments into the system except O2 on PG. Especially, 
the system of O2 on CG is antiferromagnetic (AFM). 



II. COMPUTATIONAL METHODS 

Spin-polarized density-functional theory (DFT) calcu- 
lations are performed using the Perdew-Burke-Ernzerhof 
(PBE)^^ generalized gradient approximation (GGA) for 
the exchange-correlation potential. A supercell of 6 x 6 
graphene including 72 C atoms with a doped atom sub- 
stituting a C atom and a single O2 molecule adsorbed 
onto it is constructed. With this model, the dopant con- 
centration in our calculations is ^ 1.4%. Experiments 
have shown the existence of these types of single-atom 
doped materials^^. Besides, a method to synthesize the 
transition metals doped graphene has been proposed^^. 
In the direction normal to the surface, 29.5 A length in 
the supercell is sufficient to minimize the interaction be- 
tween graphene layers. Ultrasoft pseudopotentials^^ and 
a plane-wave basis set-up to a kinetic energy cutoff of 25 
Ry for the wave function and of 200 Ry for the charge 
density are chosen in all simulations (for transition met- 
als Cr and Mn, a cutoff of 35 Ry for the wave func- 
tion and 300 Ry for the charge density). The Brillouin 
zone is sampled using a 3 x 3 x 1 Monkhorst-Pack^^ grid 
and Methfessel-Paxton^^ smearing of 0.01 Ry. A denser 
11 X 11 X 1 Monkhorst-Pack grid and the tetrahedron 
method^*" are used for the calculation of density of states 
(DOS) and partial DOS (PDOS). Atomic positions are 
optimized until the maximum force on any atom is less 
than 0.001 a.u. All calculations are performed using the 
Quantum-ESPRESSO package^^. We carefully test this 
supercell for the convergence of energy, magnetic proper- 
ties, comparing a larger supercell as in Ref. 21, and there 
is few difference found. 

Here, we would like to point out that the usage of 
GGA, the consequent neglect of van der Waals interac- 
tions, leads to an incorrect description of physisorption, 
but this is of little concern for us since we are interested 
in chemically bound molecules. In fact, in physisorbed 
systems, LDA results look "better" than GGA, but they 
are just "wrong in a different way" . For chemically bound 
systems, GGA is usually a better choice. 

In this calculation, only the triplet O2 molecule is con- 
sidered, which is the ground state of O2. The adsorption 
energy is defined by Ea = Etot - Edg - Eo, where Etot is 
the total energy of the doped graphene with a bound O2 
molecule, E^g is the energy of doped graphene and Eq is 
the energy of isolated O2 molecule. In order to minimize 




(a) Mn-doped graphene (b) O2 on MG 



FIG. 1: (Color online) The most stable configurations of MG 
with an elevation h of 1.555 A (a); and of O2 on MG with h 
of 1.652 A (b). Yellow: Carbon, Cambridge blue: Mn, Red: 
O, respectively. 



systematic errors, the same supercells and /c-point grids 
are used for all calculations. 



HI. RESULTS AND DISCUSSION 

The bond length of O2 molecule in ground state is 
1.2350 A (experimental value is 1.207 A) and the mag- 
netic moment is 2 11 b in our calculations. O2 molecule 
is placed at the top of the doping atoms, and slightly 
inclined to the graphene plane at the beginning. After 
relaxation, O2 molecule is found to be physisorbed on 
the top of B and N atoms, while chemisorbed on the top 
of Al, Si, P, Cr and Mn atoms. After relaxation, O2 
is parallel to graphene plane, which is the most stable 
configuration according to the test of our calculations. 
The configurations of MG and O2 adsorption on MG are 
shown in Fig.l, where the Mn atom protrudes out of the 
graphene plane. All configurations here are very simi- 
lar to Fig. 1, except for BG and NG with all atoms in 
the same plane and physisorbed O2. In fact, a metastable 
configuration with all atoms in one plane can be obtained 
if the dopant is set to be in the same plane of graphene, 
but it is very easy to transform to the stable one. Mean- 
while, the results of O2 adsorption on all doped graphene 
are shown in Table. I, which will be discussed detailedly 
below. 



A. adsorption of oxygen on CG and MG 

We first discuss the adsorption of O2 on CG and MG. 
The atomic and magnetic structures of CG and MG have 
been discussed in Ref. 21. As shown in Fig. 1 and Ta- 
ble. I, the elevation of Cr and Mn atoms is 1.6487 and 
1.5555 A; the Cr-C and Mn-C bond lengths are 1.8559 
and 1.8317 A; the magnetic moments are 2.00/iB and 
3.00/iB, respectively. These are in good agreement with 
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TABLE I: Summary of atomic structures of the compounds, including adsorption energy Ea (eV), the shortest bond length of 
X-C dx-c (A), shorter bond length of X-0 (ix-o(A), bond length of 0-0 (io-o(A), the elevation of the dopant atoms above 
the graphene plane h (A) (negative value means the dopants protrude out of the plane in opposite direction to O2 molecule), 
magnetic moments of the system Mb (i^b), and Lowdin charge^^ transfers from doped graphene to O2 molecule c (e~). 







Cr 


Mn 


Al 


Si 


P 


B 


N 


doped 


h 


1.6487 


1.5555 


1.7584 


1.4579 


1.4591 


0.0000 


0.0000 


graphene 


dx-c 


1.8559 


1.8317 


1.8528 


1.7611 


1.7689 


1.4794 


1.4079 




Mb 


2.0000 


3.0000 


0.0000 


0.0000 


1.0500 


0.0000 


0.0000 


O2 


Ea 


-2.6098 


-2.0918 


-1.5589 


-1.3132 


-1.0359 


-0.0232 


-0.1228 


adsorption 


h 


1.6188 


1.6520 


1.9026 


1.6607 


1.4840 


-0.2457 


-0.0541 




dx-c 


1.8159 


1.8464 


1.9284 


1.8284 


1.7501 


1.4810 


1.4066 




dx-o 


1.7935 


1.8641 


1.8770 


1.7109 


1.6275 


3.5099 


3.3196 




do-o 


1.4088 


1.3993 


1.4000 


1.5103 


1.5584 


1.2354 


1.2577 




Mb 


AFM 


1.4500 


1.0000 


0.3300 


0.0000 


1.9900 


1.8000 




c 


0.2166 


0.1865 


0.4367 


0.9184 


0.8181 


-0.0863 


0.0235 



the results using a larger supercell in Ref. 21. The un- 
paired electrons in transition metals (Cr and Mn) induce 
magnetic moments, and the direction of spin polarization 
of Cr or Mn is opposite to that of the nearest carbon 
atoms, that is, Ca, Cb and Cc in Fig. la. 

From Table. I, we can learn the configurations are 
changed dramatically after O2 adsorption. The molecu- 
lar oxygen is chemisorbed on CG and MG with large ad- 
sorption energies of -2.6098 and -2.0918 eV, respectively, 
indicating very stable adsorption of O2. The bonds of 
Cr-0 and Mn-0 are formed with short lengths of 1.7935 
and 1.8641 A; The elevation of Cr above graphene plane 
is shortened to 1.6188 A, while the elevation of Mn is 
extended to 1.6520 A. The 0-0 bond extends to 1.4088 
and 1.3993 A for CG and MG. Interestingly, CG become 
AFM by the adsorption of O2, i.e., with zero total mag- 
netic moment but 0.34/iB absolute magnetic moment. 
We calculate the magnetic moment of every atom in the 
system, and find the magnetic moment of Cr atom (- 
0.1435) is opposite to other atoms such as O (0.041 1/i^, 
0.0421/i^), the three nearest C atoms Ca (0.0133/1^), 
Cb (0.0169/iB) and Cc (0.0169/1^). At the same time, 
the absolute magnetic moment of Cr and the sum of mag- 
netic moment of O, C^, Cb and Cc is very close. These 
two properties make the system AFM, and the magneti- 
zation very localized, as shown in Fig. 2a. For MG, O2 
adsorption reduces the magnetic moment of the system 
to 1.45 /iB, which is caused by the decrease of unpaired 
electrons due to the hybridization between Mn and O 
atoms. Similar to O2 on CG, the magnetic moment is 
located around the Mn atom, as shown in Fig. 2b, which 
may be a Kondo system. 

In order to understand the electronic properties of the 
system, DOS and PDOS of O2-CG and O2-MG are cal- 
culated, as shown in Fig. 3 and 4, respectively. It is well 
known that intrinsic graphene is semimetal. In Fig. 3, it 
is obvious that CG is a metal. After adsorption of O2, 
the O, C and Cr hybridize strongly. There is one DOS 
peak for spin down electrons appearing at the Fermi level, 
which comes from the hybridization of Cr atom and O2 
molecule. One peak disappears below the position of -1.5 



eV, which mainly belongs to the contribution of 3d elec- 
trons in Cr atom. Lowdin charge analysis shows there is 
0.2166 e~ transfer from CG to O2 molecule, that is to say, 
O2 is an acceptor here. Therefore, O2 can actually p-dope 
the host CG, which is consistent with the effect of O2 on 
organic semiconductors^^. Similarly, O2 can p-dope MG 
through accepting 0.1865 e~ . From Fig. 4, we can learn 
that the change of electronic properties of MG induced 
by O2 adsorption is mainly around the Fermi level, where 
the DOS peak of spin-down electrons of MG before ad- 
sorption above the Fermi level disappears, which is con- 
tributed by the 3d electrons in Mn atom hybridizing with 
2p electrons in O atom. Meantime, one peak for spin- up 
electrons appears at the Fermi level, which is made up of 
hybridized Mnd and Op orbitals according to the PDOS 
analysis in Fig. 4. Similar to the analysis of spin den- 
sity (Fig. 2), the magnetization is also localized and then 
O2 adsorbed on Cr or Mn may exhibit as Kondo impu- 
rities as well as Cr or Mn doped into graphene^ ^ and Ni 
impurity in a Au nanowire^^. 




(a) O2 on CG (b) O2 on MG 



FIG. 2: (Color online) Spin density in O2 adsorption on doped 
graphene. (a) in O2-CG chemisorption system with isovalues 
of ±0.01 a.u.; (b) O2-MG chemisorption system with isovalues 
of dzO.OOl a.u.. Blue color represents spin-down states, red 
color spin-up states. 
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FIG. 3: (Color online) Total and partial DOS for CG and 
O2-CG chemisorption system, (a) Total DOS of CG before 
adsorption of O2; (b) Total DOS of CG after adsorption of O2; 
(c) PDOS of atoms in O2-CG chemisorbed system. PDOS of 
Op is the average of two O atoms, and Cp is the average of all 
C atoms in the system. The arrows denote the spin-down (^) 
states and spin-up (t) states. The zero energy is the Fermi 
level. 
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FIG. 4: (Color online) Total and partial DOS for MG and 
O2-MG chemisorption system, (a) Total DOS of MG before 
adsorption of O2; (b) Total DOS of MG after adsorption of 
O2; (c) PDOS of atoms in O2-MG chemisorbed system. The 
arrows denote the spin-down (X) states and spin-up (t) states. 
The zero energy is the Fermi level. 



B. adsorption of oxygen on AG, SG and PG 

Al, Si and P elements contain 3p orbital, and they 
have 3, 4 and 5 valence electrons respectively. Discussion 
about their doping into graphene and nanotubes has been 
reported recently, which show that they are potential re- 
sources for detecting toxic molecules and modulating the 
electronic structures of graphene. 

From Table. I, we can learn that Al, Si and P atoms 
introduce a local curvature of graphene, with elevation 
of 1.7584, 1.4579 and 1.4591 A, respectively. Meanwhile, 
the bond lengths of Al-C, Si-C and P-C is 1.8528, 1.7611 
and 1.7689 A, respectively, which basically decrease with 
the increase of electrons of the elements. In fact, for 
atoms which have the same configurations of orbit als, 
the covalent bond length decreases with the increase of 
group number. This indicates that the extension of elec- 
tron states plays a very important role in the structures 
of doped graphene. Specially, P dopant introduces spin- 
polarization into the PG with magnetic moment of 1.05 
/i^, while no magnetization exists for AG and SG. It is 
noticed that a metastable configuration of PG with all 
atoms retaining in one plane can be observed, in which 
there is no spin polarization induced. It seems that the 
local curvature has a very important effect for the mag- 
netic property. O2 molecule can be chemisorbed on the 
atoms of Al, Si and P in graphene, with large adsorp- 
tion energies of -1.5589, -1.3132 and -1.0359 eV, respec- 
tively. Meantime, the atomic structures change much. 
For AG, Al atom has an elevation of 1.9026 A, and C-Al 
bond extends to 1.9284 A, Al-0 to 1.8770 A and 0-0 
to 1.40 A. For SG, Si atom protrudes out of graphene 
plane with elevation of 1.6607 A. C-Si, Si-0 and 0-0 
bonds also elongate, which are 1.8284, 1.7109 and 1.5103 
A, respectively. On the contrary, chemisorption of O2 
on PG shorten the P-C bond to 1.6286 A from 1.7691 
A before adsorption, and the bond 0-0 is broken com- 
pletely with the long distance of 1.5598 A. The P atom 
almost retains the place above graphene with 1.4796 A 
and does not protrude outward more. It is interesting 
that chemisorption of O2 introduces spin polarization for 
AG and SG with magnetic moments of 1.00 and 0.33 /i^, 
but the magnetization in PG vanishes after O2 adsorp- 
tion. For AG and SG, the spin-polarization is mainly 
located in the O atoms, which is also very localized as 
shown in Fig. 5. Therefore, O2 adsorption may introduce 
Kondo effect into AG and SG systems. 

The electronic properties are also investigated by the 
analysis of DOS (PDOS) and Lowdin charges. For AG 
after adsorption, the analysis of Lowdin charges shows 
that O2 molecule accepts about 0.44 e~ from Al atom, 
and 0.72 e~ transfers from Al to the three nearest C 
atoms (0.24 per atom). The DOS of AG before and after 
adsorption of O2 is shown in Fig. 6. Al has 3 valence 
electrons, and therefore can be hole-doping for graphene, 
as shown in Fig. 6a, where the minimum of DOS shifts 
above the Fermi level. Chemisorption of O2 introduces 
spin polarization, which creates the difference of the DOS 
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Top view 




(a) O2 on AG 



(b) O2 on SG 



FIG. 5: (Color online) Spin density in O2 adsorption on doped 
graphene. (a) in O2-AG chemisorption system; (b) O2-SG 
chemisorption system. Blue color represents spin-down states, 
red color spin-up states. The isovalues are ±0.002 a.u. 
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FIG. 6: (Color online) Total and partial DOS for AG and 
O2-AG chemisorption system, (a) Total DOS of AG before 
adsorption of O2; (b) Total DOS of AG after adsorption of 
O2; (c) PDOS of atoms in chemisorbed system. The arrows 
denote the spin-down (l) states and spin-up (t) states. The 
zero energy is the Fermi level. 
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FIG. 7: (Color online) Total and partial DOS for SG and 
O2-SG chemisorption system, (a) Total DOS of SG before 
adsorption of O2; (b) DOS of SG with 3% concentration; (c) 
Total DOS of SG after adsorption of O2; (d) PDOS of atoms 
in chemisorbed system. The arrows denote the spin-down (^) 
states and spin-up (t) states. The zero energy is the Fermi 
level. 
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mainly around the Fermi level, as shown in Fig. 6b. The 
peak at Fermi level for spin down electrons is contributed 
by the Op orbital according to the PDOS analysis, shown 
in Fig. 6c. 

For chemisorption of O2 on SG, O2 molecule accepts 
about 0.92 e~ from Si atom, and about 0.23 e~ transfers 
from Si atom to every nearest C atom. Si has 4 elec- 
trons, the same as C atom. Therefore, the doping of Si 
should not shift the Fermi level of graphene, as shown 
in Fig. 7a. It is worth pointing out that the electronic 
structures are dependent on the dopant concentration, as 
shown in Fig. 7b, and the electronic band of SG system 
is open with a gap of 0.1 eV. Nevertheless, it does not af- 
fect the configuration of chemisorption of O2 on SG. O2 



FIG. 8: (Color onhne) Total and partial DOS for PG and 
O2-PG chemisorption system, (a) Total DOS of PG before 
adsorption of O2; (b) Total DOS of PG after adsorption of 
O2; (c) PDOS of atoms in chemisorbed system. The arrows 
denote the spin-down (^) states and spin-up (t) states. The 
zero energy is the Fermi level. 



adsorption introduces magnetization with 0.33 /iB- As 
shown in Fig. 7c-d, the difference of the DOS between 
spin up and spin down electrons is mainly caused by Op 
orbital. The DOS peak of spin up electrons at the Fermi 
level is caused by the p^ orbital in C atom, while the peak 
of spin down by Op. Moreover, the hybridization of Sip, 
Op and Cp orbit als also happens here. 
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For O2 on PG, about 0.82 e~ transfers from P atom 
to two O atoms, and about 0.24 e~ from P atom to ev- 
ery nearest C atom. Before O2 binding, PG is magnetic 
with 1.05 \iB' Here, charge analysis shows us that ev- 
ery nearest C atom accepts about 0.18 e~ from P atoms, 
and about 4.3 e~ remains in P atoms. That is to say, 
there is about one electron unpaired, inducing about 1 
\iB magnetic moment. After adsorption of O2, the hy- 
bridization between P and O, even inluding C, employs 
the unpaired electron, and then the spin polarization van- 
ishes. According to the DOS of PG before adsorption of 
O2, as shown in Fig. 8a, the difference of DOS between 
spin up and spin down electrons mainly locates around 
the Fermi level. The binding of O2 evaporates the mag- 
netization and creates a peak at the position of -0.75 eV, 
as shown in Fig. 8b, which is caused by the hybridization 
of Pp, Op and Cp orbitals, as shown in Fig. 8c. 



C. Adsorption of oxygen on BG and NG 

B and N belong to the second group elements in the 
periodic table, which have 2s2p configurations for elec- 
trons and are the nearest elements for C. The extension 
of electrons of B, C and N is similar. Therefore, there is 
no local distortion when B and N are doped in graphene. 
The length of B-C is 1.4794 A, larger than that of N- 
C with 1.4079 A. B and N atoms do not have elevation 
above graphene, and there is also no spin polarization. 
O2 is physisorbed on BG and NG, which is the same as 
O2 on B-doped graphite^^. The physisorption does not 
change the structures of BG and NG significantly, and 
the magnetic moments of O 2 is also almost retained, as 
shown in Fig. 9a. The spin density of O2-NG is very 
similar. The spin density of O2-BG system concentrates 
in the O2 molecule, which is almost the same as the spin 
density of isolated O2 molecule, shown in Fig. 9b. The 
O2 molecule is far away from graphene with a distance 
larger than 3.5 A and a relatively small adsorption en- 
ergy. Furthermore, B and N atoms is slightly pushed be- 
low graphene plane by O2, indicating the repulsed force 
between them. Since B and Al, N and P belong to the 
same group elements and have very similar valence elec- 
tron configurations, and their behaviors caused by O2 
molecule are so different, it can be induced that one im- 
portant reason for this phenomenon is the local curvature 
in AG and PG. This is also in agreement with the con- 
clusion of local curvature enhancing chemical reactivity 
of doped graphene^^. About this local curvature, as in 
Table. I, only the B and N dopants are not elevated, and 
the B-C and N-C bond lengths are close to C-C bond in 
graphene. By contrary, other dopants are elevated, and 
the X-C bonds are much longer than C-C bonds. There- 
fore, it can be concluded that dopant X is elevated only 
when the length of X-C is appreciably larger than that 
of C-C. 

Although O2 can not bind stably on BG and NG, it can 
still affect the electronic properties at low temperature at 




(a) O2 on BG 



FIG. 9: (Color online) Spin density in O2 adsorption in (a) 
O2-BG physisorption system, and (b) isolated O2. The iso- 
values are dbO.005 a.u. 
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FIG. 10: (Color online) (a) Total DOS of BG before adsorp- 
tion of O2; (b) Total DOS of BG after adsorption of O2; (a) 
Total DOS of NG before adsorption of O2; (d) Total DOS of 
NG after adsorption of O2; Red solid line represents spin- up 
states and blue dashed line spin-down states in (b) and (d). 
The zero energy is the Fermi level. 



least. The dopant of Boron can hole-dope graphene and 
Nitrogen electron-dope graphene, as shown in Fig. 10a, c, 
where the Dirac point shifts above (BG) or below (NG) 
the Fermi level. After adsorption of O2, the spin polar- 
ization is introduced and the DOS is changed. One DOS 
peak around Fermi level appears for spin down electrons 
of BG and NG physisorbed O2 molecule, caused by Op 
orbital; but the DOS of spin up electrons is very simi- 
lar to that of BG before adsorption. This means that 
the change in DOS is concentrated on the O2 molecule. 
However, this effect should not exist at room tempera- 
ture because of the very weak adsorption and the long 
distance between O2 and graphene. Furthermore, the 
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charge transfer between O2 and BG or NG is very small 
compared with the chemisorption, indicating the small 
change for the conductance of the system. It should be 
noticed that O2 acts as a donor for BG through donating 
~0.08 e~ to BG system. Although the charge transfer 
of the physisorbed system using the figure of the single 
electron Kohn-Sham equation is under debate, the trend 
of this small charge transfer is reliable. 

IV. CONCLUSION 

In conclusion, our simulation studies reveal that most 
of doped graphene are sensitive to O2 molecule. The lo- 
cal curvature caused by doped atoms is very important 
to the reactivity for gas molecules. BG and NG are in- 
ert to molecular oxygen, while graphene with Al, Si, P, 
Cr and Mn doping allows O2 to form stable chemisorp- 
tion states, which affects the magnetic, electronic and 
atomic properties of graphene. O2 acts as an acceptor 
in all chemisorbed configurations, and introduces local- 
ized spin polarization except on PG. In particular, O2 on 
CG is AFM, and O2 can retain its magnetic moments on 
B- and N-doped graphene. Significantly, O2 adsorption 



may introduce Kondo effect into CG, MG, AG and SG 
systems. Most of systems with O2 on doped graphene 
behave as metallic materials. Therefore, a combination 
of foreign atom doping followed by exposure to air may 
be an effective way to tune the electronic and magnetic 
properties of semimetal and unpolarized graphene. How- 
ever, the potential usage of doped graphene as gas sensors 
should be dependent on the sensitivity to O2 molecule in 
air. When chemisorption of O2 on doped graphene hap- 
pens, the sensitivity to other gases should be affected 
significantly, which may prevent it from being effective 
gas sensors. 
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